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ABSTRACT 


An  account  is  given  of  an  experiment  which  toe  place 
in  Alberta,  Canada,  during  the  summer  of  1963 .  The  primary- 
purpose  of  the  experiment  was  to  look  for  correlation  between 
ionospheric  phase  height  fluctuations  and  long  period  mag¬ 
netic  micropulsations  with  the  thought  that  this  correlation 
might  be  interpreted  either  in  terms  of  hydromagnetlc  waves 
or  as  drift  motions  produced  by  E  region  currents.  Atmo¬ 
spheric  pressure, fluctuations  were  also  recorded  for  possible 
correlation  studies. 

The  experimental  equipment,  techniques,  and  parameters 
are  all  described  and  consideration  is  given  to  those  aspects 
of  the  experiment  which  might  merit  modification  in  any  sim¬ 
ilar  experiment  in  the  futurp. 

The  data  are  compared  in  three  ways: 

1.  Phase  height  —  micropulsations 

2.  Phase  height  —  atmospheric  pressure 

3.  Ordinary  wave  --  Extraordinary  wave  phase  height 

Although  positive  correlation  was  not  observed,  it  is 
shown  how  the  lack  of  correlation  could  be  attributed  to  cer¬ 
tain  definite  causes  which  can  be  eliminated  in  any  future 
experiment.  • 
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I.  INTRODUCTION 


An  interesting  aspect  of  the  study  of  the  ionosphere  is 
the  consideration  of  wave  propagation  and  structure.  Hines 
(i960)  considers  wave  structure  observed  in  rocket  trails 
and  interprets  these  observations  in  terms  of  internal  at¬ 
mospheric  gravity  waves.  Using  radio  techniques,  Chan  (19b. ) 
observed  traveling  ionospheric  disturbances  (TID's)  and  on 
several  occasions  noticed  strong  correlation  between  TID's 
and  fluctuations  in  the  earth's  magnetic  field.  Rishbeth 
and  Garriott^  (1964)  consider  Chan's  observations  and  inter¬ 
pret  the  correlation  both  in  terms  of  the  "dynamo-motor 
theory  and  in  terms  of  ionospheric  hydromagnetlc  waves. 

Any  wave,  whatever  its  character  and  origin,  should  cause 
motion  of  the  plasma  as  it  propagates  through  the  ionosphere. 
It  would  appear  then,  that  a  method  for  sensing  small,  li^al- 
lzed  oscillations  in  the  ionosphere  would  be  an  effective 
tool  for  the  study  of  a  variety  of  wave  phenomena  associated 
with  the  ionosphere. 

A  simple  but  highly  sensitive  cw  phase  height  sounder  has 
been  developed,  and  the  experiment  herein  described  was  set 
up  around  this  device  in  an  attempt  to  study  the  hypotheses 
noted  above.  Briefly,  the  technique  is  as  follows:  l)  Ver¬ 
tical  motion  of  the  ionosphere  is  measured  by  means  of  the 
phase  sounder.  2)  Magnetic  field  fluctuations  are  measured 
along  three  components.  3)  Atmospheric  pressure  fluctuations 
are  measured.  Following  the  work  of  Hines,  we  look  for  cor¬ 
relation  between  l)  and  3).  Following  the  work  of  Rishbeth 
and  Garriott,  we  look  for  correlation  between  l)  and  2). 

It  is  instructive  to  consider  in  greater  detail  the  two 
types  of  waves  which  have  been  postulated,  i.e.,  hydromag¬ 
netlc  waves  and  gravity  waves. 
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A .  HYDROMAGNETIC  WAVES 


A  great  deal  is  known  about  hydromagnetlc  waves  from 
studies  in  plasma  physics.  In  the  paper  by  Rishbeth  and 
Garriott,  only  the  simplest  form  of  hydromagnetlc  wave  is 
considered,  namely  an  Alfven  wave  propagating  along  a  mag¬ 
netic  field  line.  On  the  basis  of  this  and  certain  other 
assumptions,  an  expression  is  derived  which  relates  (vertical) 
plasma  velocity  to  magnetic  field  fluctuations.  The  expres¬ 
sion  is 

VZH  =  cot  1  ('“Ain  V"iN^  (^ 

in  which  is  the  vertical  velocity  of  the  plasma,  AB 

is  the  amplitude  of  the  magnetic  field  fluctuation,  I  is 
the  dip  angle,  vL  is  Ionic  collision  frequency,  \iQ  is 
the  permeability  of  free  space,  mL  is  ionic  mass,  and  N 
lr  electron  number  density.  The  significance  of  Eq.  (l), 
f->r  purposes  of  this  experiment,  lies  in  its  prediction  of 
proportionality  between  vertical  plasma  velocities  and  the 
•  amplitudes  if  magnetic  field  fluctuations .  In  Section  III, 
this  relationship  is  explored  quantitatively,  and  expressions 
are  derived  predicting  the  relative  amplitudes  of  the  para¬ 
meters  to  be  measured  in  this  experiment. 


D.  GRAVITY  WAVES 

The  existence  of  atmospheric  gravity  waves  may  be  postu¬ 
lated  by  considering  the  conditions  which  govern  atmospheric 
oscillations.  These  conditions  can  be  expressed  in  terms  of 
three  equations:  the  equations  of  motion,  of  adiabatic  state, 
and  of  continuity.  The  simultaneous  solutions  of  these  equa¬ 
tions  predict  the  propagation  of  waves.  Hines  has  solved 
these  equations  and  shown  that  where  gravity  is  neglected. 
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the  solutions  reduce  to  the  ordinary  solutions  of  sound  pro¬ 
pagation,  and  are  therefore  called  "acoustic"  waves.  Under 
these  conditions  the  equation  relating  wave  frequency  and 
speed  Is  non  dispersive,  i.e.,  =  C‘  Kc ,  in  which  K  Is 

the  associated  wave  number  and  C  Is  the  speed  of  sound. 

In  the  presence  of  gravity,  however,  the  solutions  give 
rise  to  permissible  values  for  u>  which  do  not  represent  a 
continuous  spectrum,  but  rather  fall  into  two  distinct  re¬ 
gions.  Only  values  of  u>  falling  within  these  regions  can 
represent  propagating  waves.  Frequencies  falling  within 
these  two  regions  are  called  and  and  are  specified 

as : 


Moreover,  since  Y,  which  is  the  usual  ratio  of  specific 
heats,  is  always  less  than  2,  then  u-a  >  always,  is 

associated  with  acoustic  waves  and  u;  is  associated  with 

&  II 

another  sequence  of  waves  which  have  been  labeled  as  in¬ 
ternal  gravity  waves." 

In  principle,  values  for  tt-  can  be  arbitrarily  small. 

© 

Indeed,  observations  have  indicated  the  existence  of  waves 
with  peri  ds  of  the  order  of  30  minutes  for  which  the  inter¬ 
pretation  in  terms  of  gravity  waves  poses  interesting  experi¬ 
mental  possibilities. 


-  3  - 
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II.  DESCRIPTION  OF  EXPERIMENT 


A.  IONOSPHERIC  PHASE  SOUNDING 

The  phase  soundex*  measures  vertical  motion  of  the  iono¬ 
sphere  by  measuring  the  modulating  effect  this  motion  has  on 
a  signal  received  from  a  source  of  px’eclslon  frequency. 
Bi‘iefly,  the  operation  is  as  follows:  A  cw  signal  is 
transmitted  from  a  source  which  is  driven  by  a  frequency 

O 

standard  (accuracy  one  part  in  10  per  day).  The  resulting 
sky  wave  is  reflected  from  the  ionosphere  and  received  at  a 
recording  site  a  short  distance  away  (Figure  l).  The  receiv¬ 
ed  signal  is  compared  with  a  reference  signal  controlled  by 
another  frequency  standard  and  the  instantaneous  phase  dlf- 
fer'ence  is  recorded.  In  the  absence  of  ionospheric  motion, 
the  summation  of  the  received  signal  and  the  reference  signal 
will  x-esult  In  a  constant  frequency  difference.  If  the  re¬ 
ference  signal  is  adjusted  such  that  tne  difference  frequency 
Is  -ex*o,  then  the  phase  difference  will  be  unchanging  with 
t  lmc  . 

Undex*  these  clx-cumstances,  the  Doppler  effect  of  vertical 
motion  of  the  plasma  will  result  in  a  phase  modulation  of  the 
received  signal  which  may  be  observed  as  a  variation  in  the 
signal  rcpi’esenting  the  phase  difference. 

A  detailed  explanation  of  the  system  operation  is  divided 
into  four  parts:  Transmitter  Operation,  Ionospheric  Effects, 
Receiver  Operation  and  Signal  Processing,  and  Phase-Height 
Pai-ameters .  These  will  now  be  considered  in  oi’der. 

1 .  Tx-ansmltter  Operation 

A  self-explanatory  block  diagram  of  the  transmitter 
system  is  shown  in  Figure  2.  The  four  available  frequencies 
were  2.0,  2.2,  4.0,  and  4.4  Me.  Two  antenna  systems  were  em¬ 
ployed  with  half  wavelength  center  frequencies  of  2.1  Me  and 

-  4  - 
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transmitter 

BUILDING 


PLAN  VIEW  OF  EXPERIMENT  SITE 


4.2  Me.  The  antenna  impedance  at  2.0  Me  turned  out  to  be 
lower  than  expected,  and  this  resulted  in  excessive  loading 
of  the  transmitter.  Consequently,  2.0  Me  was  used  very  little. 


FIG.  2.  PHASE  SOUNDER  TRANSMITTER  SYSTEM 
.  Ionospheric  Effects 


A  sketch  of  Idealized  ray  paths  is  shown  in  Figure  3. 
As  shown  in  the  sketch,  the  received  signal  may  be  considered 
as  a  summation  of  a  sky  wave  S^  and  a  ground  wave  Sg. 

Both  Sj^  and  S^  are  of  the  form 

P 

S  =  A  cos(wt  - J  3dp)  (4) 
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SITE  SITE 


FIG.  3.  PHASE  SOUNDER  PROPAGATION  PATHS 


as  seen  by  the  receiving  antennae.  Here  dp  Is  the  lncre- 

2n 

mental  propagation  path  length  and  3  =  -y  where  X  Is  the 
wavelength.  For  the  ground  wave,  the  expression  becomes 

S,  =  A2  cos ( itt  -  SgD)  (5) 


in  which  3  =  - —  where  X  is  the  free  space  wavelength 

0  xo 

(constant).  For  the  sky  w ave,  3  Is  a  function  both  of 

space  and  time.  However,  we  can  define  the  refractive  Index 
of  a  medium  as  ^  =  X ^/X  so  that  0  =  ^3,.  Lhe  sky  wave 
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signal  is  then  given  by 


fp 

Si  =  Ai  cos(uot  -  3  I  ^ 


The  total  signal  S  seen  by  the  receiving  antennas  is  there¬ 
fore  given  by 

S  =  Sl  +  S2  (7) 


Under  typical  conditions,  the  2.0  Me  and  2.2  Me  sky 
wave  signals  tended  to  be  absorbed  during  daylight  hours. 

Under  these  conditions,  k1  «  A2  and  S  «=  S,  .  This  pro¬ 
vided  a  convenient  means  for  calibrating  the  frequency  stand¬ 
ards.  Except  for  this,  however,  the  presence  of  a  ground  wave 
was  a  nuisance.  In  general,  useful  data  was  obtained  only 
when  A^  «  Ai  so  that  S  ^ 

We  can  greatly  simplify  the  problem  and  yet  not  lose 
very  much  accuracy  if  we  assume  vertical  incidence.  This  is 
approximately  true  as  shown  by  Figure  3.  On  this  basic,  we 
can  define  a  quantity  "phase  height"  by  the  expression 

1 ' 

H1  =  2  u  dp  =  u  dh  where  H  =  P/2  and  is  considered  to  be 

Jo  Jo 


the  true  height  of  reflection  of  the  wave.  We  may  now  express 


S1  ar'  po 

S1  =  A1  cos ( xt  -  2So  I  udh )  I 

U  O 

The  frequency  of  is  given  by  ( a-  -  ^2SoJ  ond 

since  x  is  essentially  constant,  changes  in  frequency  are 

A  1 


due  entirely  to  time  variations  in  ( 


Unfortunately, 


the  relationship  is  not  easily  expressed  without  making  fur¬ 
ther  simplifying  assumptions.  Specifically.,  assumptions  must 
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be  made  with  respect  to  p.  as  a  function  both  of  time  and  of 
height.  The  assumption  we  make  which  greatly  simplifies  the 
problem,  but  is  reasonable  from  a  physical  point  of  view,  is 
as  follows:  Ionization  is  not  a  continuous  function  of 
height,  but  exists  only  in  a  layer  with  a  well  defined  lower 
boundary.  A  change  in  phase  is  caused  by  a  vertical  motion 
of  the  entire  layer.  The  shape  of  the  layer  does  not  vary 
with  time. 


If  we  designate  the  bottom  of  the  layer  as  H^,  we 


can  say 


udh  = 


udh  + 


From  the  above  assumption 


OH 


and 


fH 

I  udh 

jHb 

nHb  f 

,  I  udh 

Jo  Jo 


(9) 


dh  =  H, 


udh  is  constant;  therefore,  =  H  minus  a  con¬ 


stant.  We  can  now  express  as 


•  Sx  =  AICOs(u;t  +  -  2£qH  ) 


(10) 


where  0^  is  a  constant. 

3  •  hecc-lvrr  Operation  and  Signal  Processing 

Opera tl  n  of  the  receiving  system  may  best  be  under¬ 
stood  by  referring  to  Figure  4.  The  signal  S  =  S1  =  A.^ 

c  os ( u  t  -  25  ,11  +  is  received  by  the  master  receiver  as 

shown.  The  reference  signal  =  A,  cos(u)t  +  0^)  Is  -In¬ 
jected  into  the  slave  phase  receiver  by  virtue  of  extracting 
an  appropriate  harmonic  from  a  precision  100  kc  signal.  At 
this  point,  the  phase  difference  between  and  is 

simply  A0  $  -  +  2eQH.  and  are  individually 

-  9  - 
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FIG.  4.  TYPICAL  OF  TWO  RECEIVING  SYSTEMS 


heterodyned  down  to  a  nominal  80  kc  In  two  stages.  This  is 
accomplished  with  common  oscillators  having  isolated  outputs 
so  that  the  80  kc  signals  retain  the  precision  phase  ulff  r- 
ence  ^  -  9±  +  2SQ  H.  At  this  point,  the  signals  are 

individually  (frequency)  divided  by  four  to  nominal  values 
of  20  kc  which  is  within  the  operating  frequency  of  the  Acton 
phasemeter.  The  phasemeter  output  is  a  voltage  proportional 


-  f 


y.  e  h 


Slnce 


^  -  9, 


is  arbitrary  and 


can  be  suppressed,  we  see  that  the  time  vai*ying  portion  of 
the  phasemeter  is  a  direct  analog  of  H(t).  As  noted  in  the 
figure,  division  by  eight  is  also  available,  allowing  two 
basic  levels  of  system  sensitivity.  Details  of  sensitivities 
and  scale  factors  are  considered  in  another  section  and  will 
not  be  further  discussed  here.  Note  that  in  the  presence  of 
a  strong  ground  wave,  S.S2=A2  co«(.t  +  P2)  and  the  phase¬ 
meter  output  provides  a  true  indication  or  the  difference  in 
frequency  between  the  frequency  standards.  The  presence  of 
a  dally  ground  wave  permitted  continuous  checks  on  stability 

of  the  standards. 

An  additional  feature  of  the  system  was  the  Doppler 
or  offset  frequency  comparison.  The  offset  frequency  stand¬ 
ard  war  tuned  to  provide  a  reference  signal  (Gi()  several 
cyclcc-per-sccond  away  from  S^.  This  signal  is  received 
and  processed  in  a  manner  similar  to  that  described  above. 

The  1-f  form  of  can  be  expressed  as  A^cosU^t  +  jfi) 

while  the  1-f  form  of  S]L  is  of  the  form  A^osUt  -  29QH). 
The  sum  of  these  two  signals  is  detected  and  fed  into  a  slow 
speed  tape  recorder.  The  recorder  input  has  the  form 
A  C0stU4  -  U.)t  +  230ii  +  Y)  =  A  COsUa-t  +  °)  •  iG  Vhc  (Prc" 
cislon)  frequency  difference  resulting  from  off-setting  one 
standard  and  Is  the  frequency  of  a  precision  carrier  which  is 
phase  modulated  by  variations  in  H.  Am  was  made  as  small 
as  possible,  consistent  with  the  low  frequency  response  limit¬ 
ations  of  the  recorder. 
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4 .  Phase-Height  Parameters 


Phase  height  (i.e.,  phase  path  length)  was  measured 
as  a  function  of  time.  There  were  eight  possibilities  for 
received  signals;  ordinary  and  extraordinary  waves  on  each 
of  the  four  frequencies:  2.0,  2.2,  4.0,  and  4.4  Me.  Two 
phase  height  parameters  were  received  and  recorded  simultan¬ 
eously  in  various  pair  combinations.  The  2. 0/2. 2  Me  and  the 
4. 0/4. 4  Me  combinations  were  to  provide  redundancy  in  the 
event  of  broadcast  interference.  By  this  technique,  it  was 
intended  that  a  "low"  frequency  and  a  "high"  frequency  be 
kept  available  for  the  experiment .  The  equipment  did  not 
function  well  at  2.0  Me  and  occasional  broadcast  interference 
was  experienced  at  4.4  Me,  so  it  turned  out  that  almost  all 
of  the  records  obtained  were  combinations  of  2.2  Me  and  4.0 
Me  signals.  In  addition,  the  daylight  signal  on  2.2  Me  con¬ 
tained  a  strong  component  of  ground  wave.  Data  acquisition 
tended  to  adhere  to  the  following  pattern:  l)  During  day¬ 
light  hours  and  into  the  early  evening,  recorded  parameters 
were  usually  4.0  Me  ordinary  and  extraordinary  phase  height. 
2)  When  the  critical  frequency  of  the  ionosphere  passed  be¬ 
low  4.0  Me,  the  equipment  was  switched  over  to  2.2  Me  ordin¬ 
ary  and  extraordinary.  3)  When  this  signal  reverted  back  to 
a  ground  wave  sometime  after  sunrise,  the  equipment  was 
switched  back  to  4.0  Me  again.  4)  On  come  occasions,  var¬ 
ious  combinations  of  2.2  Me  versus  4.0  Me  signals  were  re¬ 
corded  simultaneously. 


B.  MAGNETIC  MEASUREMENTS 

The  magnetic  measurements  were  made  entirely  with  equip¬ 
ment  provided  by  Pacific  Naval  Laboratories  at  Victoria,  B.  C. 
The  parameters  required  for  the  experiment  were:  l)  Total 
magnetic  field,  and  2)  Two  components  of  magnetic  micro- 
pulsations,  true  north  (X)  and  true  East  (Y).  These  will  be 
considered  in  order. 

SEL-b4-06l 
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1.  Total  Magnetic  Field 


Total  magnetic  field  was  measured  by  a  rubidium  vapor 
magnetometer.  Functioning  of  the  equipment  was  erratic  and 
the  output  signal  was  very  noisy  to  the  extent  that  signal 
sensitivity  was  limited  to  about  three  gammas.  This  defect 
was  not  prohibitive  since  the  total  field  was  largely  used 
as  a  qualitative  check  on  the  micropulsation  signals.  The 
magnetometer  signal,  along  with  the  X  and  Y  components 
of  micropulsation  signals  provided  sufficient  information  to 
determine  amplitudes  and  polarizations  of  magnetic  mlcropuls- 
atlons . 

2 .  Magnetic  Micropulsations 

The  magnetic  micropulsation  parameters  were  true  (geo¬ 
graphical)  North  (X)  and  true  (geographical)  East  (Y).  The 
transducers  employed  were  of  the  Inductor  type,  and  as  such, 
did  not  have  a  uniform  response  over  the  frequencies  of  inter¬ 
est.  The  response  of  the  equipment  peaked  sharply  at  about 
0.5  cps,  and  fell  off  at  about  20  db  per  decade  at  lower  fre¬ 
quencies.  The  effect  of  this  non-uniformity  was  that  system 
gain  had  to  b<?  adjusted  for  the  stronger  signals  appearing  at 
0.5  cps  with  the  result  that  low  frequency  signals  were  ex¬ 
tremely  weak.  The  effect  of  this  problem  on  interpretation 
of  results  is  discussed  more  fully  in  the  next  section. 

C.  ATMOSPHERIC  PRESSURE  MEASUREMENTS 

A  microbarograph  was  set  up  to  record  atmospheric  pressu  ^e 
waves  in  the  region  of  one  millibar  amplitude  and  2-20  minute 
periods.  The  actual  range  scales  employed  were  0.02  psi  and 
0.064  psi,  and  the  low  frequency  limit  of  the  equipment  was 
estimated  to  be  (20  minutes)-1.  Associated  with  the  signal 
output  of  the  equipment  was  a  drift,  probably  thermal  in  nature, 
with  a  peak  to  peak  amplitude  of  about  0.1  psi  and  a  period  of 
about  24  hours. 
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FIG.  5.  INTERIOR  OF  RECORDING  SITE 
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FIG.  6.  TRANSMITTER  SITE 


FIG.  7.  INTERIOR  OF  TRANSMITTER  BUILDING  SHOWING  FREQUENCY 
STANDARD,  FREQUENCY  SYNTHESIZER  AND  TRANSMITTER. 
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FIG.  8.  RECEIVING  SITE 


III.  INTERPRETATION  OF  RESULTS 


The  experiment  yielded  a  considerable  amount  of  data. 

It  was  necessary  therefore  to  limit  the  scope  of  the  analysis 
to  those  areas  of  the  data  most  likely  to  yield  interesting 
and  useful  results.  For  this  reason,  the  interpretation  of 
results  was  limited  to  three  major  areas;  l)  Phase  height 
and  micropulsations,  2)  Phase  height  and  atmospheric  pres¬ 
sure,  and  3)  Ordinary  wave  and  Extraordinary  wave  phase 
height.  The  first  of  these  represents  the  primary  justifica¬ 
tion  for  the  experiment;  consequently,  most  of  the  effort  has 
been  expended  in  this  area.  The  three  areas  are  considered 
separately  in  the  following  sections. 


A.  PHASE  HEIGHT  AND  MICROPULSATIONS 

A  plot  of  ordinary -wave  and  extraordinary -wave  phase 
height  and  the  peaks  of  the  'X'  component  (N  -  S)  of  micro¬ 
pulsations  is  shown  in  Figure  9  for  a  three-hour  period  on 
July  16,  1963.  Approximately  sixty  percent  of  the  micropuls¬ 
ation  peaks  appeared  to  correlate  with  phase  height  data. 

This  correlation  occured  early  in  the  experiment  and  was  very 
encouraging.  Unfortunately,  this  degree  of  correlation  was 
not  attained  at  any  other  time  during  the  course  of  the  ex¬ 
periment  . 

In  attempting  to  establish  correlation  between  phase 
height  fluctuations  and  magnetic  micropulsations,  it  is  in¬ 
structive  to  consider  the  relative  amplitudes  of  signals  one 
should  expect  to  observe  on  the  records.  The  derivation  of 
an  appropriate  prediction  involves  essentially  two  steps; 
l)  prediction  of  the  relative  characteristics  of  the  physical 
quantities  themselves  and  2)  description  of  the  measuring 
equipment  and  its  effect  on  the  measured  quantities.  These 
steps  will  now  be  considered  in  order. 
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1.  Physical  Quantities 

The  most  conclusive  observation  of  these  phenomena  ap¬ 
pears  to  be  the  experiment  described  in  the  report  by  Chan. 

As  noted  In  the  report,  correlation  was  observed  between  mag¬ 
netic  mlcropulsatlons  and  frequency  excursions  about  a  standard 
frequency,  which,  under  certain  assumptions,  can  be  interpreted 
in  terms  of  vertical  motion  of  a  region  of  the  ionosphere. 

In  their  paper,  Rishbeth  and  Garriott  derive  three  in¬ 
dependent  expressions  for  vertical  velocity.  Of  these,  VZR 
is  shown  to  follow  directly  from  the  assumptions  and  approxima¬ 
tions  made  by  Chan  (l.e.,  mirror  reflection,  two-hop  propaga¬ 
tion,  etc.)  and  is  linearly  related  to  the  observed  frequency 
excursions.  Two  other  expressions,  VZM  and  VZR,  are  dev¬ 
eloped  from  theory  and  each  shows  a  linear  relationship  between 
vertical  velocity  and  magnetic  micropulsations.  In  Table  2  of 
their  paper,  these  quantities  are  tabulated  and  compared  with 
some  of  Chan's  observations.  The  table  is  repeated  here  for 
easy  reference.  An  additional  parameter,  Vzr/aB,  has  also 
been  included. 


Table  2 


Table 

1  Ref. 

Approx . 
Local 
Time 

hr. 

Period 
of  Osc. 

T. 

min. 

Freq . 
Devn. 
Af 
c/s 

Horlz . 
Mag.  Devn. 
AB  (Fr) 

V 

v 

VZR 

m/s 

VZM 

v 

VZH 

VZR 

AB 

E 

11 

.56 

0.3 

1.4 

5 

0.6 

1.8 

3.56 

G 

12 

0.7 

0.4 

0.4 

5 

0.2 

0.6 

12.5 

A 

13 

2.1 

0.6 

2 

8 

0.9 

1.0 

4.0 

B 

15 

3.2 

0.8 

3.2 

11 

2.0 

1.5 

3.44 

F 

17 

4.7 

1.1 

2.4 

8 

5 

1.0 

3.33 

H 

18 

6 

0.8 

1.5 

6 

4.5 

0.8 

4.0 

C 

19 

1.3 

0.5 

0.5 

6 

1.0 

0.6 

32.0 

D 

21 

3.7 

0.7 

2 

10 

12 

1.8 

5.0 

J 

22 

15 

0.6 

1.5 

8 

9 

0.6 

5.33 
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Examination  of  the  values  listed  for  Vzr/AB  (m/s  psr  y) 
reveals  that  seven  of  the  nine  values  are  reasonably  self 
consistent  with  a  mean  of  about  4  m/s  per  y,  more  or  less 
independent  of  the  perturbation  period.  With  the  further 
assumption  that  Chan's  micropulsation  equipment  had  a  "flat" 
response  over  the  range  of  frequencies  Involved,  we  may  con¬ 
clude  that  the  ratio  of  amplitudes  of  the  physical  quantities 
themselves  does  not  vary  with  perturbation  period.  For  the 
purposes  of  this  exercise,  the  value  4m/s  per  y  will  be 
considered  representative. 

Based  on  the  foregoing,  we  can  say 


AP  =  2(helght  change)  =  2  I  VZR  DT 


J. 


(11) 


where  At  is  the  time  interval  of  Interest,  and  subscript 
"a"  refers  to  actual  values. 

If  we  assume  (for  simplicity)  sinusoidal  variations,  and 
let  ABa  =  ABa  cos  ait,  then  V^R  =  4aB3  cos  (tut  +  0) 

f  8aB  (12) 

and  APa  =  2  J  4ABa  cos  (tut  +  0)dt  =  — —  sln(tut  +  0) 

Disregarding  phase,  the  ratio  of  peak  values  (Fg)  is  given  by 

AP  o  .  \ 

F  =  — —  =  —  meters/ganuna  (13) 

a  AB„  tu 

a 


Moreover,  if  we  let  tu  =  2™/!,  then 

AP  h*!1  /  I  \ 

F  =  =  —  meters/gamma  (1^ ) 

a  ABa 

Subject  to  the  various  assumptions  made,  Fa  is  the  actual 
quantity  of  interest. 

2 .  Measuring  Equipment 

There  were  essentially  three  pieces  of  equipment  pertinent 
to  this  portion  of  the  experiment;  the  phase  sounder,  the  micro- 
pulsation  gear,  and  the  Sanborn  recorder. 
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a)  The  phase  sounder  output  was  a  signal  proportional 
in  amplitude  to  the  path  length  measured  in  units  of  the 
wavelength  of  the  operating  frequency.  As  part  of  the  de¬ 
tailed  mechanics  of  the  system,  the  signal  was  routed  through 
a  unit  which  performed  a  "by  4"  or  "by  8"  division  at  all 
times.  The  phase  sounder-recorder  combination  was  set  up  to 
have  the  following  scale  factor  (S  ): 

s' 


Sp  =  £§  divisions/meter 


(lb) 


where  X  =  wavelength  of  operating  frequency 
N  =  signal  division  factor  (4  or  8) 


so  that 


4o( AP  ) 

AP  =  (S  )( AP  )  =  — a 
r  '  p/v  a'  \N 


(lb) 


where  subscript  "r"  refers  to  recorded  values.  For  the  usual 
conditions  of  4  Me  operating  frequency  and  N  =  4,  we  have 
X  =  75  meters,  and 

APr  =  S( AP^ )  =  (  AP.  )  =  y- r  Sanborn  record  divisions  ( 17 ) 

l)  The  scale  factor  for  the  Sanborn  recorder-micropulsa¬ 
tion  gear  combination  (Sm)  was  determined  from  a  frequency 
response  calibration  of  the  equipment.  The  inverse  scale 

factor  Gm  =  -^-  is  plotted  in  Fig.  10.  Then 
in  kj 

m 

AH 

AH  =  (S  )(  ABJ  =  — p—  Sanborn  record  divisions  (l8) 

r  '  m  a  u 

m 


Although  we  cannot  conveniently  describe  G  analytically, 

AP  m 

I* 

we  can  express  the  ratio  F^  =  — g—  . 
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PERIOD  IN  SECONDS 


FIG.  10.  MICROPULSATION  SCALE  FACTOR 


AP 


F  = 
r 


40(APa) 


m 


40G 


m 


AP„ 


40G  F0 
m  a 


lbOG  T 
m 


A  B_ 


\N 


AB, 


XN 


AB„ 


XN 


nXN 
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The  function  F^  describes  the  ratio  of  peak  values  of  cor¬ 
related  signals  which  one  might  expect  to  observe  on  the  re¬ 
cords,  if  such  correlated  signals  are  present  . 

Four  typical  cases  are  plotted  in  Fig.  11. 

(1)  f  =  4  Me  (x  =  75)  and  N  =  4 

(2)  f  =  4  Me  and  N  =  8 

(3)  f  =  2  Me  (X  =  150)  and  N  =  4 

(4)  f  =  2  Me  and  N  =  8 
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11.  CORRELATION  AMPLITUDE  RATIO  FROM  EXPERIMENTAL  DATA 
B.  PHASE  HEIGHT  VERSUS  ATMOSPHERIC  PRESSURE 

Ac  noted  earlier,  theory  predicts  the  existence  of  at¬ 
mospheric  pressure  waves  of  extremely  low  frequency.  As 
p  )inted  out  ty  Hines,  the  postulation  of  such  waves  provides 
an  adequate  interpretation  of  various  observed  phenomena  such 
as  TID's  and  meteor  trail  distortions.  The  gravity  waves  can 
propagate  in  any  direction  with  both  transverse  and  longitud¬ 
inal  components  of  perturbation.  The  vertical  components  of 
such  perturbations  would  result  in  vertical  oscillations  of 
the  ionospheric  plasma  which  in  principle  could  be  detected 
by  a  vertical  incidence  sounder.  The  pressure  oscillations 
associated  with  such  waves  could  be  detected  through  the  use 
of  a  microbarograph.  Unfortunately,  typical  vertical  phase 
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velocities  (tens  of  meters  per  second)  suggest  that  correla¬ 
tion  of  atmospheric  disturbance  effects  might  be  prevented 
by  the  great  time  delay  involved.  In  addition,  pressure  os¬ 
cillations  observed  at  a  single  location  provide  no  informa¬ 
tion  with  regard  to  phase  velocity,  wavelength,  or  direction 
of  wave  normal. 

The  outstanding  characteristic  of  any  wave  which  should 
be  evident  in  any  of  its  manifestations  is  that  of  frequency 
of  oscillation.  Correlation  could  exist  in  the  following 
manner;  a  propagating  gravity  wave  should  cause  plasma  motion 
at  ionospheric  heights  and  pressure  variations  at  ground  level. 

It  would  be  extremely  difficult  to  predict  with  any  degree 
of  accuracy  the  relative  amplitudes  of  the  two  parameters  one 
might  expect  to  measure.  It  is  interesting  to  note,  in  this 
regard,  that  a  pressure  measurement  detects  a  wave  through 
its  compressional  characteristics  which  are  associated  with 
longitudinal  motion  of  the  medium.  Meteor  trail  observations, 
however,  indicate  that  gravity  waves  at  ionospheric  heights 
are  characterized  primarily  by  transverse  motions  of  the 
medium. 

An  additional  factor  to  be  considered  is  that  the  origin 
of  gravity  waves  is  not  well  understood.  Without  knowing 
where  a  wave  originated,  the  time  delay  between  two  spatially 
separate  measurements  cannot  be  predicted. 

In  Fig.  12  is  shown  a  typical  plot  of  atmospheric  pressure 
variations  and  4.0  Me  Ordinary  wave  phase  height  variations, 
both  as  functions  of  time.  No  correlation  was  detected  from 
this  plot,  nor  was  any  correlation  detected  from  any  of  the 
data . 
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TIME  (Z) 


l’’IG .  12.  TYPICAL  COMPARISON  OP’  PHASE  HEIGHT  AND  ATMOSPHERIC 
PRESSURE  VARIATIONS 

C.  ORDINARY  WAVE  VS  EXTRAORDINARY  WAVE  PHASE  HEIGHT 

hrough  the  use  of  crossed  dipole  antennas  and  hybrid 
transformer  units,  the  phase  sounder  receiver  was  able  to  re¬ 
ceive  separately  and  simultaneously  both  the  Ordinary  (0)  and 
Extraordinary  (X)  components  of  the  reflected  wave.  Accord¬ 
ing  to  the  magnetoionic  theory,  the  ray  paths  of  0  and  X 
waves  diverge  upon  entering  the  ionosphere  so  that  the  re¬ 
flected  components  of  these  separate  waves  represent  reflec¬ 
tion  from  different  regions  of  the  ionosphere.  It  has  been 
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shown  by  Budden^  (1961)  and  elsewhere  that  the  ray  path  of  a 
vertically  incident  wave  lies  in  a  plane  defined  by  the  mag¬ 
netic  meridian.  In  the  case  of  vertical  incidence  which  is 
appropriate  to  this  experiment,  the  ray  paths  are  approximate¬ 
ly  as  shown  in  Fig.  13  in  which  the  reasonable  assumption  is 
made  that  collision  frequencies  are  much  smaller  than  the  wave 
frequency.  As  seen  in  the  sketch,  the  points  of  reflection 
of  the  two  waves  differ  both  in  latitude  and  elevation. 


FIG.  13.  ORDINARY  AND  EXTRAORDINARY  RAY  PATHS 
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It  is  instructive  to  consider  the  latitudinal  separation 
between  the  reflection  heights  of  the  Ordinary  and  the  Extra¬ 
ordinary  wave.  If  we  assume  vertical  incidence  and  no  colli¬ 
sions,  we  can  make  use  of  equation  13.44  in  Budden  as  follows: 


H  =  * 


sin  0  cos  0  Yc( 

7 -  - 


4Yt2(1 


X)2  + 


(20) 


in  which  0  is  the  angle  between  the  ray  path  and  the  verti¬ 
cal,  1  is  a  horizontal  coordinate  in  the  magnetic  meridian, 
and  z  is  a  vertical  coordinate. 

If  we  let  dl  equal  the  incremental  lateral  deviation  of 

a  ray,  then  we  can  say  dl  =  (-|i)dz  and 


for  the  ordinary  wave  (2l) 


for  the  extraordinary  (22) 
wave 


where  Z  and  Z  are  the  heights  of  reflection  for  the  re- 
cpective  waves. 

Substituting  ( 20 )  into  (2l)  and  (22),  we  obtain 


L  = 


v 


Z  sin  0  cos  0  Y‘  (l  -  n  ) 

o  _ _ _ o  ' 

o  ^4  Yl2(1  -  X )2  +  YT4j  2 


dz 


L  = 


Zv  sin  9  cos  9  Y‘"(l  -  n^  ) 

x)J  +  yt4J  2 


Ic  sin  a  cos  * 
o  ^4  Yl2(1  - 


dz 


(23) 


(24) 
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We  will  use  a  Chapman  model  ionosphere,  so  for  convenience, 
we  make  a  vertical  translation  of  our  coordinate  system  so 
that  vertical  distance  is  measured  with  the  peak  of  the  layer 
considered  as  the  origin.  Because  in  the  Chapman  model  the 
electron  concentration  is  negligible  at  distances  greater  than 
three  scale  heights  below  the  peak,  the  integration  lower  limit 
is  set  at  minus  three  scale  heights.  A  reasonable  value  for 
daytime  scale  height  is  80  km,  so  our  lower  limit  becomes  240 
km.  We  then  have 


-240 


Z  sin  9  cos  0  Y2(l  -  n  c  ) 
°  Q 

4  Yl2(1  -  X)2  +  Y TUJ  5 


r 


dz 


=  P 


Z„  sin  0  cos  0  Y‘‘(l  -  ny  ) 


2  .  v  M  * 

T 


(• 


2H0  l"  yL  <1  '  X>"  +  Y' 


dz 


(25) 


(26) 


where  z  is  a  coordinate  directed  upward  from  the  peak  of  the 
Chapman  layer.  The  various  functions  appearing  in  the  inte¬ 
grals  are  standard  magnetolonlc  parameters  and  are  defined  be¬ 
low  as  by  Ratcliffe5  (1959): 


Y 


ui  m 


(27) 


Y^  =  Y  cos  0 


(28) 


YT  =  Y  sin  0 


(29) 


X 


Ne 


2 


E  m 
o 


2 


a' 


(30) 
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From  the  Appleton-Hartree  equations,  we  have 


1 


1 


o 


(31) 


The  Chapman  layer  is  defined  by 


N 


(33) 


The  value  of  Nmax  is  determined  from  the  following  con¬ 
siderations:  (l)  Reflection  of  the  ordinary  wave  occurs 

where  X  1,  that  is  where  the  plasma  frequency  of  the  med¬ 
ium  becomes  equal  to  the  v/ave  frequency.  (2)  Associated  with 
the  layer  is  some  maximum  frequency  at  which  the  ordinary  wave 
can  be  reflected  at  vertical  incidence.  If  we  designate  this 
frequency  as  fmx  and  define  «max  -  2"(  fmax),  then  Nmax 
is  determined  directly  from  equation  (30)  as: 


X 


1  = 


N  e 
max 


e  m(u)  Y 
o  max ' 


hence 


N 

max 


e  m(  x’  )L 
o  max' 


(34) 


(3b) 
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From  the  condition  that  the  Ordinary  wave  reflects  at 
X  =  1  and  the  Extraordinary  wave  reflects  at  X  =  1  -  Y 
(see  Ratcliff e)  we  can  determine  ZQ  and  Z x,  the  upper 
limits  of  our  two  integrals.  To  determine  ZQ,  we  set  X  =  1 
and  z  =  Zq  and  combine  equations  (30),  (33)  and  (33), 

whence ; 


from  which  is  obtained 

+  exp(  -  -pp)  =  1  -  2 log  — ■ — — ~2  (36) 

n  (  O)  ) 

v  max' 


Similarly,  for  Z x  we  have  X  =  1  -  Y,  z  =  Zx,  and 


which  leads  to 


Equations  (36)  and  (37)  each  have  two  solutions;  that  is, 

for  any  — —  <  1,  each  equation  has  one  positive  and  one 
“■'max 

Z 

negative  solution  for  .  The  physical  reason  for  this  is 

that  for  any  value  of  N  (N  <  ^max)  In  the  layer,  there  is 

one  point  above  the  peak  and  one  point  below  the  peak  where 
such  a  value  of  N  can  be  found.  Since  we  are  concerned  only 
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with  the  region  below  the  peak,  only  negative  values  of 
are  considered. 


The  geography  of  the  experiment  determined  the  values 
9  =  13°  and  m.qH0  =  0.601  x  10  4  webers/m^.  Using  a  radio 
frequency  of  4.0  Me  and  a  variety  of  representative  values 
for  fmax'  solutions  to  equations  (25),  (26),  (36)  and 

(37)  have  been  obtained  through  the  use  of  a  high  speed  dig¬ 
ital  computer.  Integrals  {25)  and  (26)  have  been  solved 
in  iterative  form,  yielding  a  complete  ray  path  profile  for 
each  combination  of  parameters,  although  only  the  values  for 
Z  ,  Z  ,  L  ,  and  L  are  of  immediate  interest. 

vJ  A  O  A 

These  results  have  been  tabulated  in  Table  2  where  values 
of  the  quantities  have  been  found  for  various  values  of  f 

max 

and  a  fixed  wave  frequency  f  =  4.0  Me. 


Table  2 


The  situation  is  a  realistic  one:  variations  in  N 

max 

give  rise  to  a  proportionate  variation  of  f  ,  according 
to  equation  (35).  It  can  be  seen  that  for  transmitting  fre¬ 
quencies  very  close  to  the  critical  frequency,  the  quantity 
(L  +  L  )  is  quite  sensitive  to  changes  in  the  critical  fre- 
quency.  An  accurate  determination  of  the  parameter  (L^  +  L^) 
would  require  the  use  of  an  ionosound  or  other  means  to  provide 
a  continuous  measurement  of  f  . 
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A  close  inspection  of  Pig.  14  reveals  that  the  two  phase 
height  curves  are  generally  similar  but  that  much  of  the  fine 
structure  is  not  coincident  in  time.  This  suggests  traveling 
waves  and  leads  to  the  conclusion  that  this  technique  can  be 
used  as  a  means  of  studying  traveling  disturbances  in  the 
ionosphere.  Davies^  (1962)  has  commented  on  this  technique 
and  noted  the  precautions  which  should  be  observed  in  the  in¬ 
terpretation  of  results. 


IV.  CONCLUSIONS 


As  noted  earlier,  the  experiment  had  three  aspects,  which 
were  essentially  Independent.  The  comparison  of  phase  height 
with  micropulsations  was  the  primary  justification  for  the  ex¬ 
periment.  The  other  two  aspects,  comparison  of  phase  height 
with  atmospheric  pressure  and  comparison  of  Ordinary  wave  and 
Extraordinary  wave  phase  height,  were  pursued  insofar  as  they 
did  not  Interfere  with  the  primary  consideration.  Neverthe¬ 
less,  conclusions  can  be  drawn  for  each  part  of  the  experiment, 
and  these  will  now  be  considered  in  order. 


A.  COMPARISON  OF  PHASE  HEIGHT  AND  MICROPULSATIONS 


The  purpose  of  this  portion  of  the^  experiment  was  to  look 
for  correlation  between  ionospheric  phase  height  fluctuations 
and  long  period  magnetic  micropulsations.  Such  correlation 
was  observed  in  an  experiment  by  Chan,  and  an  interpretation 
of  the  correlation  has  been  suggested  by  Rishbeth  and  Garriott 
as  noted  earlier.  Inasmuch  as  observed  correlation  has  been 
rare,  the  experiment  herein  described  was  intended  to  provide 
experimental  evidence  of  this  interesting  phenomena. 

The  observations  were  inconclusive.  Positive  correlation 
was  not  observed,  yet  it  is  impossible  to  conclude  that  such 
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FIG.  14.  ORDINARY  AND  EXTRAORDINARY  PHASE  HEIGHT, 
JULY  1963 
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correlated  physical  phenomena  were  necessarily  absent.  In 
examining  the  various  factors  which  Influenced  the  experi¬ 
mental  results,  two  of  these  factors  have  emerged  as  having 
a  predominantly  deteriorating  effect  on  the  experiment. 

These  factors  are  as  follows:  l)  The  response  of  the  mag¬ 
netic  micropulsation  equipment  was  not  matched  to  the  require¬ 
ments  of  the  equipment.  The  effect  of  this  mismatch  is  ex¬ 
plained  in  Section  III-C,  so  it  will  not  be  repeated  here. 

It  might  be  said  that  the  narrow  bandwidth  associated  with 
the  equipment  used  in  this  experiment  appears  to  be  character¬ 
istic  of  the  coil  or  inductance  type  of  instrument.  In  con¬ 
trast,  the  micropulsation  gear  used  in  the  Chan  experiment 
had  a  uniform  response  to  all  frequencies  of  interest. 

2)  The  radio  equipment  did  not  function  properly  at  2.0  or 
2.2  Me.  The  transmitter  antennas  were  forty  feet  from  the 
ground,  and  the  results  of  this  low  height  were  excessive 
loading  of  the  transmitter  and  excessive  ground  wave  emis¬ 
sions.  The  ground  wave  interfered  with  the  night-time  sky 
wave  with  the  result  that  very  little  useful  data  could  be 
obtained  on  either  2.0  or  2.2  Me. 


B.  COMPARISON  OF  PHASE  HEIGHT  AND  ATMOSPHERIC  PRESSURE 

Of  the  three  aspects  of  the  experiment,  this  appears  to 
be  least  promising.  The  biggest  drawback  appears  to  be  that 
a  given  pressure  fluctuation  at  ground  level  which  is  thought 
to  be  associated  with  a  gravity  wave  provides  no  information 
with  regard  to  direction  or  velocity  of  the  gravity  wave. 
Since  the  phase  sounder  is  exploring  only  a  small  region  of 
the  ionosphere,  the  chance  of  finding  correlation  appears 
remote.  This  portion  of  the  experiment  was  performed  only 
because  it  was  so  easily  accomplished j  it  involved  only 
setting  up  and  recording  the  output  of  a  microbarograph. 
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C.  COMPARISON  OF  ORDINARY -WAVE  AND  EXTRAORDINARY -WAVE  PHASE 

HE-lGHT 

Although  this  aspect  of  the  experiment  has  not  been  dev¬ 
eloped  very  far,  the  possibilities  for  further  development 
and  experimentation  appear  very  good,  indeed.  Under  the  pre¬ 
sent  circumstances,  there  are  several  drawbacks  which  should 
be  corrected  if  this  area  of  investigation  were  to  be  pursued, 
l)  The  ionospheric  model  used  here  is  a  Chapman  layer  which 
is  reasonably  valid  only  under  conditions  of  equal  production 
and  loss.  2)  The  parameter  (Lq  +  Lx)  is  of  great  importance 
and  should  be  determined  as  accurately  as  possible.  Inspec¬ 
tion  of  Table  A  reveals  the  dependence  of  this  parameter  on 
f  .  For  this  reason,  fmav  should  be  determined  at  all 
times  when  measurements  are  being  made.  3)  A  knowledge  ox 

(L  +  L  )  and  a  measurement  of  the  time  between  similar 
'  o  x / 

occurrences  on  the  0  and  X  records  will  provide  informa¬ 
tion  about  one  velocity  component  of  a  traveling  disturbance, 
this  component  being  along  a  line  connecting  the  points  of 
reflection,  that  is  (ZQ,  LQ )  and  (Zx,  1^).  If  two  receiving 
stations  were  set  up  on  a  line  normal  to  the  magnetic  mer¬ 
idian,  such  an  arrangement  would  provide  monitoring  of  four 
points  in  the  ionosphere  instead  of  only  two.  These  four 
points  would  very  closely  define  a  plane  with  the  result  thxf 
the  velocity  of  any  traveling  disturbance  could  then  be  des¬ 
cribed  in  two  dimensions  instead  of  one. 
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